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1. Introduction 


An important goal of early vision is the computation of a representation of the 
visible surfaces in an image, in particular the determination of the orientation of 
those surfaces as defined by their local surface normals [Marr 1982, Brady 1982]. 
Many processes contribute to achieving this goal, stercopsis and structure from 
motion being currently the most studied in image understanding. Three other 
important contributing processes arc shape-from-contour, shape-from-texturc-- 
gradients, and shape-from-shading. Several psychophysical demonstrations show 
that shape-from-contour is significantly more powerful than shape-from-texture- 
gradients fClark et. al. 1956, Gruber and Clark 1956, Braunstein and Payne, 
1969], Similarly, Barrow and Tcnenbaum [1981, Figure 1.3 ff] suggest that 
shape-from-contour is a more effective clue to shape than shape-from-shading. 

In this paper we consider the computation of shape-from-contour. Figure 1 
shows a number of shapes that are typically perceived as images of surfaces 
which are oriented out of the picture plane. The corresponding (sets of) surface 
normals (up to tilt reversal) arc shown to the right. It may be supposed (see 
for example [Gregory 1973, pages 168ff]) that the slant judgements in Figure 
1 are largely determined by familiarity with regular shapes such as circles and 
squares. Figure 2 strains that hypothesis (though Gregory' proposes that we are 
familiar with the shape of puddles) and suggests that the computation is based 
on more general knowledge of shapes and surfaces. The method we propose 
is based on such general knowledge, namely a preference for symmetric, or at 
least compact, surfaces. Note that the contour docs not need to be closed in 
order to be interpreted as oriented out of the image plane. Finally, Figure 
3 shows that, in general, contours are interpreted as curved three-dimensional 
surfaces. 

We develop an extremum principle for determining three-dimensional surface 
orientation from a two-dimensional contour. Initially, we work out the 
extremum principle for the cases illustrated in Figures 1 and 2, that is, assuming 
a priori that the contour is closed and that the interpreted surface is planar. 
Later, we discuss how to extend our approach to open contours and how to 
interpret contours as curved surfaces as shown in Figure 3. 

The extremum principle maximizes a familiar measure of the compactness or 
symmetry of an oriented surface, namely the ratio of the area to the square of 
the perimeter. It is shown that this measure is at the heart of the maximum 
likelihood approach to shape-from-contour developed by Wilkin [1981] and 
Davis, Janos, and Dunn [1982]. The maximum likelihood approach has had 
some success interpreting irregularly shaped objects. The method is ineffective, 
however, when the distribution of image tangents is not random, as is the 
case, for example, when the image is a regular shape, such as an ellipse or 



Figure 1 

Two-dimensional contours that are often interpreted as planes that are oriented with 
respect to the image plane. The commonly judged slant is shown next to each shape. 



Figure 2 

Some unfamiliar shapes that are also interpreted as planes that are oriented with respect 
to the image plane. The shape on the left is from [Gregory 1973, fig. 10.9], the others 
from [Witkin 1980, p 29 and 94], 


a parallelogram. Our extremum principle interprets regular figures correctly. 
We show that the maximum likelihood method approximates the extremum 
principle for irregular figures; but that the maximum likelihood method docs 
not compute the correct slant for an ellipse. Witkin [1981, Figure 5] provides 



Figure 3 

Some shapes that are interpreted as curved three-dimensional surfaces. 


empirical evidence that the maximum likelihood method computes a good 
approximation to the perceived tilt but underestimates the slant. We prove in 
Appendix A that the maximum likelihood method consistently overestimates the 
slant of an ellipse.' A more thorough investigation of the difference between 
the Extremum Principle and the Maximum Likelihood method is needed. 

One class of figures that are readily perceived as lying in a plane other than 
the image plane are skew symmetries, which are two-dimensional linear (affine) 
transformations of real symmetries. Kanade [1981, page 424] has suggested a 
method for determining the three-dimensional orientation of skew-symmetric 
figures, under the "heuristic assumption" that such figures are interpreted as 
oriented real symmetries. We prove that our extremum principle necessarily 
interprets skew symmetries as oriented real symmetries, thus dispensing with 
the need for any^euristic assumption to that effect. Kanade shows that there 
is a one-parameter family of possible orientations of a skew-symmetric figure, 
forming a hyperbola in gradient space. He suggests that the minimum slant 
member of the one-parameter family is perceived. In the special case of a real 
symmetry, Kanade’s suggestion implies that symmetric shapes arc perceived as 
lying in the image plane, that is having zero slant. It is clear from the ellipse 
in Figure 1 that this is not correct. Our method interprets real symmetries 
correctly. 

First, we review the maximum likelihood method. In Section 3, we discuss 
several previous extremum principles and justify our choice of the compactness 







measure. In Section 4, we derive the mathematics necessary to extremize 
the compactness measure, and relate the extremum principle to the maximum 
likelihood method. In Section 5, we investigate Kanade’s work on skew 
symmetry. One approach to extending the extremum principle to interpret 
curved surfaces, such as that shown in Figure 3, is sketched in Section 6. In 
the final section, we relate this work to the psychophysical literature on slant 
estimation and image understanding work on shape from texture. 

2. The Sampling Approach 

Witkin [1981] has treated the determination of shape-from-contour as a problem 
of signal detection. Recently, Davis, Janos, and Dunn [1982] have corrected 
some of Witkin’s mathematics and proposed two efficient algorithms to compute 
the orientation of a planar surface from an image contour. Witkin's approach 
uses a geometric model of (orthographic) projection and a statistical model of 
(a) tire distribution of surfaces in space (statistics of the universe) and (b) of 
the distribution of tangents to the image contour. Wc shall adopt the geometric 
model, but dispense with the statistical model in favor of an extremum principle. 

First, the geometric model. Assume that the image plane is horizontal with 
coordinates {x,y) (see Figure 4) . To obtain a plane with slant a and tilt 
r, we rotate (x, y) by r in the image plane and then rotate the image plane 
by a about the new y axis. Assume that the coordinate frame in the plane 
(a,r) is chosen so that it projects into [x,y). (In Section 4, we describe this 
transformation more precisely, see also Davis, Janos, and Dunn [1982, p3].) 

Now suppose that a curve is drawn in the plane (o, r) and denote by /? the 
angle that the tangent makes at a typical point on the curve. Let a be the 
tangent angle in the image plane at the point corresponding to /?. Then a and 
P are related by: 


, , tan B „. 

tan(ct — t) — -. (2.1) 

coso 

We now turn to the statistical model, which consists of two assumptions 
called isotropy and independence. Isotropy reasonably supposes that all surface 
orientations are equally likely to occur in nature and that tangents to surface 
curves are equally likely in all directions. More succinctly, it is assumed that the 
quantities (/?, o, r) are randomly distributed and their joint probability density 
function ("density") D(/3,o,r) is given by [Davis, Janos, and Dunn, 1982]: 

D((3, o,t) = ~ sin o (2.2) 

7 

Wc assume that the ranges of the angles are: 0 < c < J,0 < r < 7r,0 < 
P < ir. Similarly, the density of ( o,r ) is given by : 
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Figure 4 

The geometry of orthographic projection. See text for details of notation. 




D(a, t) — — sin a 

7T 


(2.3) 


The independence assumption requires that the image tangents {a;: 1 < i < n} 
are statistically independent. That is, it is assumed that the tangent directions 
at different points on the image curve arc independent, fhis is only true if 
tire contour is highly irregularly shaped, or if the number of samples is small. 
In any case, the assumption of independence is an inherent weakness of the 
sampling approach (see for example [Witkin 1981, p. 36]). 

Witkin [1981, p. 25 - 26] shows that the joint density function D(a,a,r) is 
given by: 


D{a,a, t) 


sin a cos a 


^ cos 2 (a — r) -f cos 2 csin 2 (a — r) 
For the conditional density D{a\o,r) we find: 


(2.4) 


, , 1 coscr 

D{ol\o,t) — -. (2.5) 

n cos 2 (o; — r) -j- cos 2 osin 2 (a — r) 

Denote the sample («!, a^, ... a n ) by A (the sample is independent by 
assumption). It has conditional density 
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(T\ 




O. 


D(A\o,t) = JJ D{ai\o,r) 
1=1 


By Bayes’ formula wc obtain 


( 2 . 6 ) 


D[a r | A) = _ ^a _ 

/ / D{A\o, t)D(o, r)dodT 


(2.7) 


Observe that the numerator is independent of a and r. The sampling approach 
takes a random sample A and defines the most likely orientation of the plane 
(<r,r) to be that which extremizes D{o,t\A). Witkin [1981] quantizes a and 
r, and describes an algorithm to find the maximizing [aj,Tk). Davis, Janos, 
and Dunn [1982] develop a more efficient algorithm that first estimates a and 
r and then uses those estimates in a Newton iterative process. They provide 
evidence that their method is more accurate than Witkin’s. Curiously, however, 
they state [Davis, Janos, and Dunn 1982, p 24] that "the iterative algorithm was 
not used [in the experiments they report] because the intitial estimates (whose 
computation is trivial) are very accurate and the iterative scheme often failed 
to converge to the solution". 


3. Extremum Principles. 


Brady and Horn [1983] survey the use of extremum principles in image 
understanding. The choice of performance index or measure to be extremized, 
and the class of functions over which the extremization takes place, are justified 
by appealing to a model of the geometry or photometry of image forming and 
constraints such as smoothness. For example, the use of extremum principles 
in surface reconstruction is based upon surface consistency theorems [Grimson 
1981, 1982, 1983, Yuille 1983] and a thin plate model of visual surfaces [Brady 
and Horn 1983, Terzopoulos 1983]. 

There are several plausible measures of a curve that might be extremized in 
order to compute shape-from-contour. First, § K 2 ds, where k is the curvature 
of the contour, has been investigated as a curve of least energy for interpolating 
across gaps in plane curves [Horn 1981]. Here we seek a measure of a curve 
that is extremized when the plane containing the curve is slanted and tilted 
appropriately. Contrary to what appears to be a popular belief, given an ellipse 
in the image plane, / K 2 ds is not extremized in the plane that transforms the 
ellipse into a circle. Appendix B contains a proof of this assertion. Since 
ellipses are normally perceived as slanted circles, we reject the square curvature 
as a suitable measure. 

Another possible measure is proposed by Barrow and Tcnenbaum [1981, p89]. 
Assuming planarity (the torsion r is zero), it reduces to 
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One objection to this measure is that it involves high-order derivatives of the 
curve. This means it is overly dependent on small scale behaviour. Consider, 
for example, a curve which is circular except for a small kink. The circular 
part of the curve will contribute a tiny proportion to the integral even when 
tire plane containing tire curve is rotated. The kink, on the other hand, will 
contribute an arbitrarily large proportion and so will dominate tire integral no 
matter how small it is compared with the rest of the curve. This is clearly 
undesirable. For example, it suggests that the measure will be highly sensitive 
to noise in the position and orientation of the points forming the contour. 

A second objection to the measure proposed by Barrow and Tenenbaum is 
that it is minimized by, and hence has an intrinsic preference for, straight lines, 
for which d/c/ds zero. This means that the measure has a bias towards planes 
that correspond to the (non-general) side-on viewing position. These planes are 
perpendicular to the image plane and have slant 7r/2. 

We base our choice of measure on the following observations. 

1. Contours that are the projection of curves in planes with 
large slant are most effective for eliciting a three-dimensional 
interpretation. 

2. A curve is foreshorted by projection by the cosine of the slant 
angle in the tilt direction, and not at all in the orthogonal direction. 

We conclude that three-dimensional interpretations are most readily elicited 
for shapes that are highly elongated in one direction. Another way to express 
this idea is that tire image contour has large aspect ratio or is radially asymmetric. 
The measure we suggest will pick out the plane orientation for which the curve 
is most compact and most radially symmetric. Specifically our measure is 


M =. 

[Perimeter) 2 


(3.1) 


This is a scale invariant number characterizing the curve. For all possible curves 
it is maximized by the most compact one, a circle. By compact, we mean most 
radially symmetric. This gives the measure an upper bound of 1/47T. Its lower 
bound is clearly zero and it is achieved for a straight line. It follows that our 
measure has a built-in prejudice against sidc-on views for which the slant is 
7r/2. 

In general, given a contour, our extremum principle will choose the orientation 
in which the dcprojccted contour maximizes M. For example an ellipse is 
interpreted as a slanted circle. The tilt angle is given by die minor axis of the 
ellipse. It is also straightforward to show that a parallelogram corresponds to a 
rotated square. Appendix C discusses the interpretation of several simple shapes. 
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In particular, an ellipse is interpreted as a slanted circle, a parallelogram as a 
slanted square, and a triangle as a slanted equilateral triangle. In Section 5 we 
extend the parallelogram result to the more general ease of skewed symmetry. 

We note that the quantity M is commonly used in pattern recognition and 
industrial vision systems [Agin 1980, Pavlidis 1977, Ballard and Brown 1982] 
as a feature that measures the compactness of an object. Furthermore, we can 
show that the measure M defined in Eq. (3.1) is at tire heart of the geometric 
model in tire maximum likelihood approach. 

From Section 2, we see that the maximum likelihood approach maximizes 
the product of a number of terms of form 


/(“) 


cos a 


(3.2) 


cos 2 (cc — r) -f- cos 2 asin 2 (a — r) 

Differentiating Eq. 2.1 with respect to the arc length sj along the image curve 
and s R along the rotated curve respectively we obtain 


kj dsj 1 

K, R ds R f(a) 


(3.3) 


where k/ and k r are tire curvature at corresponding points of the image contour 
and its deprojection in the rotated plane respectively. In fact, «/ = da/dsi and 
k r = dp/ds R . There is no o or r dependence in the numerator of equation 
(3.3). We can write each term k ds as (ds ds)/(p ds) where p is the radius of 
curvature. Now observe that (pds)/(ds ds) is just a local computation of area 
divided by perimeter squared! Hence maximizing each f(a) in the maximum 
likelihood approach is equivalent to locally maximizing area over perimeter 
squared. In section (4) we will examine this connection more rigorously. 

Finally we note that the area, as well as the perimeter, can be obtained by 
an integral round tire contour. If n is the normal to the curve then it is a 
straightforward application of Stokes’ Theorem to show that 


(Area) n = - j r X dr, (3.5) 

2 4 

where (Area) is a scalar quantity, and r is a vector coordinate system in the 
plane of the figure. This formula simplifies the calculations and means the 
perimeter and the area can be computed simultaneously. 


4. Extremizing the Measure 


We now write down the measure for a curve with arbitrary orientation and then 
extremize with respect to the orientation. Let the unit normals to the image 
plane and the rotated plane be k and n respectively. The slant o of the rotated 
plane is given by the scalar product 
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cosa = k • n (4.1). 

Let T/j and L/ be the contour in tire rotated and image planes. A vector r in 
die image plane satisfies r • k — 0, and is the projection of a vector v in the 
rotated plane that satisfies v • n = 0. The projection relationship between v and 
its image r is defined by: 

r = k X (v X k) = v — (v • k)k (4.2) 

v - n x ( r x k ) _ (n • ~ (n • r)k . . 

(n-k) (n-k) 

where X denotes vector product. Now T r and F; have (vector) areas Ar and 
A/ given by 


Aj? = 

If v X dy 

(4.4) 


2 JYr 

A/ = 

1 /* 

- f r X dr 

2 Jr, 

(4.5) 


Observe that the area vectors have the same direction as the normal to the plane 
containing the area. In particular, A/? is in direction n and A/ is in direction 
k. Substituting Eq. 4.3 into Eq. 4.4 and using Eq. 4.1, we find 


It follows that 


A r = A / + 


k X (n X A/) 


(k-A/) i 
cos a 


(n-k) 


(4.6) 




M 

cos a 


(We recall that the range of a guarantees that cos a is positive. 
The perimeter lengths Pr and Pi are given by 


(4.7) 


p R = j Tr IIHI 

(4.8) 

p ‘ =1,11*11 

(4.9) 

Substituting Eq. 4.3 into Eq. 4.8 gives 



(4.10) 
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n 


In general there is no simple relationship between the perimeters analogous to 
Eq. 4.6 between the areas. Nevertheless, by Eq. 4.7, 


M = INI 

P R PrCOSG 


(4.11) 


and so our extremum principle is equivalent to extremizing cos? oPr, which 
we write as 


t 1 



(4.12) 


We extremize this with respect to the orientation n of the rotated plane, 
maintaining the constraint that n is a unit vector by a Lagrange multiplier A/4. 
This gives 


-f An = 0 

Taking scalar products of Eq (4.13) with k and n, respectively, gives 




(n • k)(n • dr)dr 
(n • k) 2 


(4.13) 



0 = (n • k)A 



(n • dr) 2 \ 2 (n • k) 2 dr 2 — (n • dr) 2 
(n • k) 2 


0 = A 

+ (n-k)jff(n.k)dr 2 + 


(n-dr) 2 A * (n • k) 2 dr 2 + (n • dr) 2 

(n-k) J (n-k ) 2 


(4.14) 


(4.15) 


We use Eq. 4.14 to remove the integral coefficient of k in Eq. 4.13, allowing 
us to express Eq. 4.13 as a sum of tire second integral and k X (k X n) times 
A. We now use Eq. 4.15 to eliminate A from this new form of Eq. 4.13. We 
recall that the unit tangent t is defined by 


dr 

d? 


where s is arc length along the contour. It follows that 


fr\ 


dr = tds. 
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Recalling that (n • k) — cos a, we find 


2 |cos 2 a + (n • t) 2 | (t • n)dr = —k X (k X n) j> jeos 2 o + (n • t) 2 | ds 

(4.16) 


where t = dr/||dr|| is the unit tangent to the image contour. 

Let the unit vectors in the x and y directions in the image plane be i and j 
and the normal to the image plane k. By definition, cos a = k • k. dlie tilt r 
is defined by 




COS T 


sin r 


n 


siner 

jn 

sin a 


The tangent vector t and the normal n can be written: 
t = cos ai + sin a\ 

n = sin a cos ri -f sin a sin rj + cos erk 


(4.17) 

(4.18) 


(4.19) 

(4.20) 


where a is the tangent angle in the image. We now form die scalar products 
of Eq. 4.16 with i and j to obtain (with appropriate cancelling) 


2 j (cos 2 o -f- sin 2 acos 2 (a — r)) ^ cos(a — r) cos ads 


COST 


(cos 2 a -f- si 


(4.21) 


sin 2 £rcos 2 (a — r))^ds 


and 


2 J> (cos 2 a -j- sin 2 crcos 2 (a — r)) = cos(a — r) sin ads 
= sin r y(cos 2 a -j- sin 2 acos 2 (a — r))^ ds 


(4.22) 


we can rewrite diese equations, after multiplying by factors of ±cosr and 
± sin r, as: 

2 y (cos 2 o -f sin 2 acos 2 (a — r)) — ^ cos 2 (a — r)ds 

= y (cos 2 a -f- sin 2 (tcos 2 (q; — r))2ds (4.23) 

2 y (cos 2 a -f- sin 2 acos 2 (a — r))~^ cos(q — r) sin(a — r)ds 


(4.24) 



Wc use Kqs. 4.23 and 4.24 in Appendix D to determine the extremi/ing 
orientation for a skew symmetry. We can rewrite these equations in the form: 


JLi 

' 1 , 

^(cos 2 a -f- sin 2 acos 2 {a 

—- r))2rfs^ = 0 

(4.25) 

da\ 

k i J 

^COS 2 o 

9 ( 

f 1 „ 

^(cos 2 a sin 2 acos 2 (a 

— r))i(isj = 0 

(4.26) 

dr\ 

l ] 

^cos2 a 


to emphasize that they correspond to extremizing with respect to a and r. 

We can implement the extremum method directly from Eqns. (4.25) and 
(4.26) by quantizing the tangent angle a;, slant Oj, and tilt t*-, replacing the 
integral by a sum. This gives good results, even though we use fixed point 
integers in our edge finder. A multi-level search speeds the algorithm by a 
factor of ten. For large slant the ratio of the greatest to least value of the 
expression is large, and the result is numerically well-conditioned. For smaller 
slants (less than about 45 degrees in die case of an ellipse) the ratio is small 
and the result poorly conditioned, so that round-off errors can be significant. 

To conclude this Section, wc show that these equations arc similar, though not 
identical, to those obtained by the maximum likelihood method in the limit as 
the number of sampled tangents tends to infinity. To see this we recall from Eq. 
2.7 that this method involves extremizing D{A\o,t) with respect to o and r. 
Since the denominator is independent of a and r, this amounts to extremizing 
D(A\a, r)D{o, t). This is the same as extremizing logD(z4|(7, r)D(a, t). Using 
2.3, 2.5 and 2.6 wc obtain: 

n 

E = n log cos a -f- log sin a — ^ log(cos 2 (a; — r) + cos 2 asin 2 (a; — r)). 

i=i 

(4.27) 

where we have ignored factors of n which will vanish on differentiation. 
Dividing E by n and taking the limit as n tends to infinity gives: 

F = log cos a J dr — j> log(cos 2 (a — r) + cos 2 asin 2 (a — r))dr. (4.28) 
Using the identity: 

cos 2 (a — r) -f cos 2 osin 2 (a — r) = cos 2 a + sin 2 <7cos 2 (a — r) (4.29) 

gives 


F = log cos a j> dr —- y log(cos 2 a + sin 2 <7cos 2 (a — r))dr. (4.30) 
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/'“'S 


r^- 


This formula is similar to Eqs 4.25 and 4.26. Thus wc expect the Extremum 
Method to give similar results to the Sampling Method when the contour is 
sufficiently irregular. However, we can show formally that the Sampling Method 
and the Extremum Method are not equivalent. In Appendix A we show that the 
Sampling Method overestimates the slant of an ellipse. The precise discrepancy 
between die methods, and its practical consequence for computing shape from 
contour is currently under investigation. 

5. Skew Symmetry 


We now consider a more general class of shapes for which the maximum 
likelihood approach is not effective. Kanade [1981, sec. 6.2] has introduced 
skew symmetries , which arc two-dimensional linear (affine) transformations of 
real symmetries. There is a bijective correspondence between skew symmetries 
and images of symmetric shapes that lie in planes oriented to the image plane. 
Kanade proposes the heuristic assumption that a skew symmetry is interpreted 
as an oriented real symmetry, and he considers the problem of computing the 
slant and tilt of the oriented plane. 

Denote the angles between the x-axis of die image and die images of the 
symmetry axis and an axis orthogonal to it (the skewed transverse axis) by a and 
ft respectively. The orthogonality of die symmetry and transverse axes enable 
one constraint on the orientation of the plane to be derived. Kanade uses 
gradient space (p, q ) [Horn, 1977, Brady, 1982] to represent surface orientations. 
He shows [Kanade 1981, p. 425] drat die heuristic assumption is equivalent to 
requiring the gradient (p, q) of the oriented plane to lie on die hyperbola 


where 


2 2 
p\ COS 



— g^sin 2 



= — cos(a — ft) 


(5.1) 



Kanade [1981, p. 426] further proposes that the vertices of the hyperbola, 
which correspond to die least slanted orientation, are chosen within this one- 
parameter family. This proposal is in accordance with a heuristic observation of 
Stevens [1980], In die special case that the skew symmetry is a real symmetry, 
that is in the case that a — ft — ^71-/2, the hyperbola reduces to a pair of 
orthogonal lines [Kanade 1981, page 426] passing through the origin. In such 
cases the slant is zero. In other words, Kanade’s proposal predicts that real 
symmetries are inevitably interpreted as lying in die image plane, and hence 






Figure 5 

A.partiLion of a skew-symmetric object into regions for the piecewise analysis presented 
in this section. 


having zero slant. Inspection of Figure 1 shows that this is not the case. A 
(symmetric) ellipse is typically perceived as a slanted circle, particularly if the 
major and minor axes do not line up with tire horizontal and vertical. 

Although Kanade’s minimum slant proposal does not seem to be correct, 
there is evidence (for example [Stevens 1980]) for Kanade’s assumption that 
skew symmetries are interpreted as real symmetries. We now show that the 
assumption can in fact be deduced from our F.xtrcmum Principle. As a corollary, 
in Appendix D we determine the slant and tilt of any given skew-symmetric 
figure; only in special cases docs it correspond to the minimum slant member 
of Kanade’s one-parameter family. 

Instead of using F.qs. 4.25 and 4.26, we prove the result from first principles 
since this enables us to use symmetry directly. First, we partition the image 
into regions as shown in Figure 5. A typical region is shown in Figure 6a and 
we skew it by an angle a to get Figure 6b. We show that the ratio of the area 
to the perimeter squared is greatest for Figure 6a with a = 0. We show this 
is also true for tire sum of the regions and take tire limit as the width of the 
regions goes to zero to complete the proof. 

From Figure 6a we calculate the area Ai ||Aj| and the tire portion Fj of 
the perimeter of the shape to be 


Ai = 2 al -f- [b — a)l 


(5.3) 
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and 

Pi = 2 s/p~+[b^a^. 

From Figure 6b we find 


and 


j4' = 2 al cos a + (6 — a)l cos a 
— A l cos a 


(5.4) 


(5.5) 


P'i = \Jl 2 cos 2 a -f- (6 — a + l sin a) 2 + \Jl 2 cos a 2 (b — a — l sin a) 2 . 

(5.6) 


We need to show that 


cos a 1 

- < —. 

p/2 p2 

1 i 1 i 

This is equivalent to showing 


(5.7) 


(2 cos a — 1)(/ 2 + (6 — a) 2 ) 

< jl 2 + (5 — a) 2 + 2 (b — a)l sin aj j/ 2 + (b — a) 2 — 2(6 — a)l sin a j . 

(5.8) 

This condition clearly holds for 7r/2 > a > n/3. Assume therefore that 
0 < a < n/3. Squaring both sides of (5.8), we see that the condition is 
equivalent to 

l A + (b - a) 4 + l\b - a) 2 (2 + —-) > 0. (5.9) 

cos^ a — cos a 

On completing the square we see this always holds provided 

4 cos a — 1 — 3cos 2 a > 0 (5.10) 

This is so provided j < cos a < 1, that is for 0 < a < |. It follows that the 
ratio of tire area to the perimeter squared is maximized for each region when 
a = 0. 

We now partition the shape into n regions as shown in Figure 5. Let the i th 
block have area A[ and perimeter P[ when the skew angle is a, and denote 
the area by A; and the perimeter by P; when a is zero. It follows from the 
above results that 




A[ — Ai cos a 


(5.11) 
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( 0 ) 



(b) 


Figure 6 

a. A basic region before skewing it. b. The result of skewing the region shown in (a) 
through angle a. 




and 


A <- A 

v ' s> _ o * 


P'. 


12 


We conclude from Eqs. (5.11) and (5.12) that 

P > Pi 2 cos a. 


Hence we obtain 


(5.12) 

(5.13) 

(5.14) 

(5.15) 

E?_iJf " EJLi p *’ (5 ' 16) 

and taking the limit as n tends to infinity shows that our extremum principle 
interprets a skew symmetry as a oriented real symmetry. 


It follows that 


Y P? > cos aY p f 

i~l i—1 

n n 

£4 

i—1 


cos a E A <- 

i—1 


EJU A < E?=i A 


2 ’ 
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6. Strategy and Sensitivity 

In Section 3, wc suggested that radially asymmetric, elongated, or non-compact, 
contours most readily elicit three-dimensional interpretations. We then proposed 
a compactness measure A/P 2 , and suggested that perceived surface orientation 
corresponds to the slant and tilt that maximize the measure. We briefly discussed 
the role of extremum principles in computer vision (Brady and Morn [1983] 
present a fuller discussion), and criticized Barrow and Tenenbaum’s measure on 
the grounds that it is highly sensitive to noise. In Sections 4 and 5 we analyzed 
the extremum principle and showed that it corresponded closely to human 
perception of oriented planar contours, and that it interprets skew symmetries 
as oriented real symmetries. 

In this Section we return to the discussion in Section 3, concentrating on two 
additional questions concerning our measure. First, we ask whether there are 
specific reasons why a visual system should adopt the strategy of extremizing 
our measure. We suggest that the extremum principle not only determines the 
orientation of the viewed curve, but provides an estimate of the stability of the 
interpretation. We relate this to the idea of general viewpoint. Second, we 
show that the compactness measure is relatively insensitive to noise and to the 
scale at which the image is sampled. To this end we consider two cases: (i) 
adding a sinusoid to the contour, (ii) assuming that the image can be modelled 
in terms of fractals [Mandelbrot, 1982]. 

Viewing position and stability 

Consider viewing a given planar curve from the hemisphere of all possible 
directions. Consider further the way the image changes when die viewpoint 
is shifted slightly. A smooth curve will hardly change with a slight change of 
viewpoint from most viewing directions. We call diese viewing directions stable 
viewpoints. Stable viewpoints can be grouped into regions whose boundaries 
correspond to viewpoints where the contour changes rapidly for a slight change 
of viewing position. These stable regions will be quite large for images of 
smooth planar curves and smooth curved surfaces. Wc suggest that image 
contours are interpretated as curves diat arc viewed from stable viewpoints. 
This is the essence of die general viewpoint constraint in computer vision. We 
are able to estimate how stable a given viewpoint is. 

One way to find stable viewing positions is to define a similarity measure for 
viewpoints and then find the extrema of tliis measure. Sutherland [personal 
communication] has proposed AjP 2 as a similarity measure, based on empirical 
studies of animal perception. If this is correct, extremizing AjP 2 corresponds 
to finding the stable viewing positions of a contour, and, as a result, die stable 
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interpretations of a contour, defined to be those which observers at genera! 
positions arc most likely to sec. This suggests drat assuming general viewing 
position corresponds to assuming interpretations around which our compactness 
measure is extremized and where it varies slowly. 

We illustrate this idea by considering an ellipse. For an ellipse the stable 
viewing position corresponds to interpreting tire ellipse as an oriented circle. 
Even if the circle is slanted by as much as thirty degrees, the value of A/P 2 
hardly changes, as our algorithm demonstrates, and so die circle is a stable 
interpretation. If the algorithm is applied to an ellipse with large eccentricity, 
corresponding to large slant, A/P 2 changes rapidly for slight changes to the 
viewing position. So if we look at an ellipse from all possible viewing angles 
it is most likely to appear as a circle and so interpreting an ellipse as a slanted 
circle is a good strategy. 

These results are preliminary, and they will be rigorously developed in a 
further paper. 

Sensitivity of the compactness measure to noise 

The second question concerns the sensitivity of the measure to noise and to 
the scale at which an image is sampled. Since our measure involves fewer 
derivatives of the contour than measures such as that proposed by Barrow 
and Tenenbaum, it should be relatively insensitive to noise. Although our 
present algorithm and the sampling mcdiod are essentially equally sensitive, it 
is possible to develop a less sensitive algorithm to implement our measure. The 
sampling method inevitably involves calculating tangents, however, and hence 
is inherently sensitive to noise and scale. For example if the image contour 
is continuous, but not differentiable, the sampling mcdiod cannot (strictly) be 
applied. Similarly if die orientation of the tangent to the contour has a large 
high frequency component, the sampling approach will be very sensitive to 
noise. 

Consider the sensitivity to noise of the extremum principle. Clearly, the 
perimeter P depends on the scale at which the image is viewed, while die area 
A is much less dependent. It follows that the ratio A/P 2 varies with scale. 
If a sinusoid is added to a smooth contour, for example, the area will remain 
approximately the same but the perimeter will change significantly. This does 
not imply, however, that the extrema of A/P 2 vary with scale. For example, it 
is easy to show that if a sinusoid is added to the contour, the extremum of the 
area over the perimeter squared is effectively unchanged. 

Similarly, suppose that the image can be approximated by a fractal [Mandelbrot, 
1982], If wc measure the image at a small fractal scale constant l, the perimeter 
will be given by 
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P(l) = Fl {i ~ d \ (6.1) 

where F is independent of l, and d > 1 is tire fractal dimension. Provided l is 
small, the value of A is essentially independent of l, since the fractal tends to a 
limit contour, with finite area, as l goes to zero. Now we extremize A/P 2 for 
all contours which can be projected into the image contour. This is equivalent 
to extremizing \og(A/P 2 ) which, using (6.1), we write as 

log(j^) = log A — 2 log F — 2(1 — d) log /. (6.2) 

We now extremize this over a and r (Eqs. 4.24 and 4.25), and note that the 
term involving l is independent of a and t and disappear. It follows that the 
extrema of A/P 2 are independent of the scale of viewing for an image that 
can be approximated by a fractal. 

7. Interpreting image contours as curved surfaces 

Figure 3 shows a number of contours that are interpreted as curved surfaces. 
In this section we discuss one method for extending our extremum principle 
to this general case. The key observation, as it was for Witkin [1981], is that 
our method can be applied locally. To do this, we assume that tire surface is 
locally planar. At the surface boundary, corresponding to the deprojection of 
the image contour, the binormal coincides with the surface nonnal. Tie idea 
is to compute a local estimate of tire surface normal by tire extremum principle 
described in the previous sections and then to use an algorithm, such as that 
developed by Terzopoulos [1983], to interpolate the surface orientation in the 
interior of the surface. The method is closely related to that proposed by Brady 
and Grimson [1981] for perceiving subjective surfaces. 

The main question concerns how to apply the extremum principle locally. 
We are currently investigating the following approach. Consider the circle of 
curvature to the space curve that is the deprojection of the contour. One way 
to define the circle of curvature is as the best fitting local circle through three 
points (or more if one assumes noisy data and makes a least squares estimate) 
on the space curve near the point in question. 'I’he circle of curvature projects 
into an ellipse. We compute die best fitting ellipse at each point on the image 
contour, and compute from it a local estimate of the surface orientation by 
finding the slant and tilt of the corresponding circle, interpreted as the circle 
of curvature. It' is easy to show that so long as the surface foreshortening is 
differentiable, this is a good estimator of the circle of curvature and hence of 
the local surface normal. 

It requires 5 parameters to define an arbitrary ellipse. Computing the 
best fitting ellipse in the general case is a complex nonlinear problem best 
approached using a numerical descent method, though several algorithms have 



been published recently for computing best fitting ellipses [Bookstein 1979, Agin 
1981, Nakagawa and Roscnfeld 1979, Sampson 1982J. If we assume that the 
normal to the space curve is not significantly foreshortened, we can compute 
the ellipse center and major axis from the curvature of the contour at the point 
in question. The slope of the contour at that point also defines the orientation 
of tire ellipse, leaving a much simpler one-parameter problem. 

We note that perceptually the strongest local cues to surface orientation 
correspond to points of high curvature. This is consistent with our method of 
locally estimating surface orientation by fitting local ellipses.. Recall that our 
compactness measure was inspired by tire observation, for example on ellipses, 
that large slant is an effective cue to surface orientation and in the case of 
an ellipse this produces points of high curvature. In fact, we can show that 
numerical conditioning of the estimator of slant increases monotonically with 
the slant. Conversely, for straight line portions of a contour, the curvature 
is zero, and the surface is locally planar. Hence surface orientation does not 
change along tine length of the straight portion. 

One issue that remains to be studied is the interface to the surface reconstruc¬ 
tion algorithm. Consider the image of a triangle shown in Figure 1. There 
are, in general, three different perceived orientations of the triangle correspond¬ 
ing to propagating the interpretation of each of the (high curvature) corners. 
Adjacent corners give inconsistent information, and so it seems necessary to use 
a labelling approach such as that proposed by Zucker, Hummel, and Roscnfeld 
[1977]. 

7. Related work 

Gibson [1950] has argued that surface orientation is directly determined by 
certain "higher order variables" in the proximal stimulus array. What exactly 
constitutes a "higher order variable” has been the subject of extensive debate. 
Nevertheless, Gibson and his followers have proposed several such, especially 
for optical flow and texture gradients. Flock [1964, Flock et. al. 1967] have 
adopted a Gibsonian perspective on the judgement of slant in texture gradients. 
He has introduced a "higher order variable" optical slant (actually "optical 
theta" in the original) as a possible basis for monocular slant perception [Flock 
1964, Eq. 5] or, at least, as a discriminant for planarity. Flock’s work implicitly 
assumes planarity [Flock 1964, p. 381], and it assumes that the tilt direction 
has been computed previously. In fact, there arc many assumptions in Flock’s 
paper (Freeman [1965, p. 502] counts 13). Moreover, several studies [Clark 
et. al. 1956, Gruber and Clark 1956, Freeman 1965, Braunstcin and Payne 
1969] have shown that the determination of surface orientation from texture 
gradients, for example using optical slant or as in lkeuchi [1983], is less effective 
than its determination from bounding contour. 
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Figure 7 

a. A box-like figure typical of those studied by Attneave and Frost [1969], b. A sheared 
box. 


Attneave [1972, see also Attneave and Frost 1969] describes an approach 
to determining surface orientation that has similarities to that developed here. 
First, he argues for a Praegnanz, theory of perception, which stresses economy 
principles in perception. This has meant different things to different researchers. 
The Gestalt psychologists explicitly noted the link between Praegnanz and 
minimization principles (for example the soap bubble) in mathematical physics. 
Attneave [1972, p. 285] suggests that such minima may be computed by "hill- 
climbing techniques". The extremum principles discussed in Section 3 of this 
paper and surveyed in [Brady and Horn 1983, Grimson 1981, 1982, 1983, and 
Terzopoulos 1983] can be considered to be more sophisticated formalizations of 
similar intuitions. 

In fact, Attneave adopts Hochberg’s formulation of tire Praegnanz theory as 
the tendency to keep differences to a minimum. In particular he considers 
three-dimensional interpretations that equalize one or more of angles, lengths 
of edges, and surface slopes in figures such as Figure 7a. The more of these 
that are in fact equalized in a particular three-dimensional interpretation of an 
image, tire more likely that interpretation is to be chosen. In fact, the extremum 
principle developed in this paper will interpret Figure 7a correctly. It will also 
determine the shear in Figure 7b. 

Second, Attneave [1972, Fig. 4] considers the judged orientation of rhombus 
figures (see Appendix D below for analysis of this case by our extremum 



principle.) Surface orientation can be judged reliably for such shapes, more 
reliably in fact than for the box figures shown in Figure 7, so long as the 
symmetries of the rhombus do not align with the horizontal and vertical. It is 
well known that the horizontal and vertical are important for shape description, 
(see for example Brady [1983] for discussion). Note that surface slant is 
consistently underestimated. 

In order to express depth constraints, Barnard [1982] and Ikeuchi [1983] 
have proposed alternative projections to orthographic projection used in the 
development in Sections 3 and 4. The methods we have developed extend 
straightforwardly to these alternative projections. 

Barnard [1982] finds surface orientation from vanishing points of families of 
parallel edges using central or perspective projection. He considers the projec¬ 
tion of an angle, and, after the fashion of Hochberg and Attneave, heuristically 
assumes equi-angularity to arrive at a three-dimensional interpretation of a 
triangle. As noted in Appendix C, our method dispenses with the need for such 
a heuristic assumption. Similar remarks apply to Barnard’s study of curvature. 

Ikeuchi [1983] proposes a method for determining surface orientation for a 
surface that is covered with uniformly repeating texture elements. He assumes 
that the surface element corresponding to the texture clement is planar. Most 
critically it is assumed that tire shape of the texture element is known a priori. 
For a particular slant and tilt, die orientation of the projection of the known 
figure changes. Ikeuchi proposes a measure that is superficially similar to the 
symmetry measure in Section 4 to determine the slant and tilt of a texture 
element. 

Olson [1974] has studied a version of Ames’ trapezoidal window illusion. 
Surface slant is judged more accurately when die stimulus is moving consistent 
with the static three-dimensional interpretation than in die purely static case. 
Similarly, Wallach, Weisz, and Adams [1956] have shown that when an ellipse 
is rotated in the image plane about an axis passing through its center and 
normal to the image it is perceived as a spinning oriented circle, like a settling 
spinning penny. The instantaneous surface orientation of die oriented circle 
can be computed by our method. Hildreth's forthcoming PhD diesis on the 
perception of motion discusses the rotating ellipse and proposed a dicoretical 
explanation of human perception of it. 
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Appendix A: The maximum likelihood method applied to an ellipse 

Equation (4.30) gives the value F of the logarithm of the Maximum Likelihood 
estimator for surface orientation: 

F — log cos a dr — j ( log(cos 2 a -f- sin 2 acos 2 (a — r))dr. (Al) 

In this Appendix we investigate the slant computed by die Maximum Likelihood 
method for the ellipse 






a 2 b 2 


= 1 . 


It is convenient to use the standard parameterization, 


(A.2) 


x = a cos# 

y = b sin#, 0 < # < 2ir (A3) 

dr = ( a 2 sin 2 # + b 2 cos 2 6)*d6 

Pending a more thorough analysis, we assume that the tilt r = ir/2 is computed 
correctly by the Maximum Likelihood method, and restrict attention to slant. 
Witkin [1981, Figure 5] and our own computational experiments suggests that 
this is reasonable; in any case, if it is not so, die inequivalence of the Maximum 
Likelihood and Extremum Principle methods would be assured. Without loss 
of generality, we assume that a > b, and denote a/b by X > 1. 

Since r = n/2, 


cos 2 a 


1 


1 + tan 2 a 

a 2 sin 2 # ’ 

a 2 sin 2 # + b 2 cos 2 # 


since tan a — dy/dx. Also, 


vr. 


sin(a — r) = sin(a-) = — cos a, 

2 


and so Eq. (A.l) reduces to 

F = log cos a dr — j> log(l — cos 2 a sin 2 a)dr. 


(AA) 


(A.5) 


{A. 6 ) 


We are interested in extrema of F, and so we consider dF/dcr, which we 
write 
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dF 
cot a — 
do 


a 2 sin 2 0 cos 2 a — 1> 2 cos 2 9 , 
-- dr, 

a 2 sin 2 0 cos 2 a -j- b 2 cos 2 0 


1 < X 


(A7) 


Suppose that the Maximum Likelihood method is extremized at a — o'. 
The slant that is necessary to interpret the ellipse as a circle is given by 
cos o — 1/X. The Maximum Likelihood method overestimates slant if and 
only if o' > cos -1 (1/X), which is if and only if cos o' < 1/X. This is true 
if dFjda 0 for all a such that 1 < Xcostr < X. Denote Xcoso by fi, 
so that 1 < fi < X. Substituting in Eq. (A.7) and changing the limits of 
integration, we find 




/i 2 sin 2 9 — cos 2 9 
H 2 sin 2 9 -f- cos 2 9 


(X 2 sin 2 9 -j- cos 2 9)*d9 


(AS) 


We now split the range of integration into four equal intervals of size 7r/4. 
With suitable changes of variables to bring the intervals of integration to [0, f], 
we find 






fj? cos 2 9 — sin 2 9 2 
fi 2 cos 2 9 + sin 2 9 
/i 2 sin 2 9 — cos 2 9 
jj? sin 2 9 + cos 2 9 


(X 2 cos 2 9 + sin 2 9 

(X 2 sin 2 9 -J- cos 2 9)*d9 


(A 9) 


Over the range of integration, cos# > sin 9. Algebraic manipulation of Eq. 
(A.9) shows it to be greater, than zero, so that it is certainly not zero. Elence 
the Maximum Likelihood method overestimates the slant. 
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Appendix B: Why square curvature is inappropriate 
In section 3 we rejected the measure 

F = j K 2 ds ( 5 . 1 ) 

because, despite popular belief to the contrary, it does not interpret an ellipse 
as a circle. In this appendix we substantiate that claim. 

Consider the ellipse defined by Eqs. (A.2) and (A.3). The curvature k of the 
ellipse is given by 




d 2 r 


= a 2 b 2 

(a 2 sin 2 9 b 2 cos 2 0) 3 

where s is the arc length. Substituting Equation (B.3) into (B.l) we have 

r2n 2l2 

F(a, b) — - de {BA) 

0 (a 2 sin 2 9 -f b 2 cos 2 6)§ 

Suppose, without loss of generality, that the ellipse is in the image plane and 
that a > b. Suppose further that the square curvature performance index (B.l) 
is correct and interprets the ellipse as a circle lying in a plane that is slanted 
with angle a to the y-axis, where 


( 5 . 2 ) 

( 5 . 3 ) 




b— a cos a (5.5) 

Equivalently, if we set b' = 6/cos a, the measure will choose o such that 
b 1 = a. Hence if we write b = \a, and consider F(a, b) to be a function of 
X, the measure should be extremized by X = 1. We will now show' that this is 
not the case. 


F(X) 


r2ir 

JO 


X 2 


aJ ° (sin 2 9 X 2 cos 2 0)t 
Differentiating this with respect to X gives 


d9 


dF 1 f 2v 2X sin 2 9 — 3X 3 cos 2 9 I(1 
— = - / - —d9 

a 0 (sin 2 9 + X 2 cos 2 6)i 


( 5 . 6 ) 


( 5 . 7 ) 



Evaluating this expression at X = 1 gives 


dF 

3X 


x=i 


i r 2n 

- (2 

a J o 


sin 2 9 — 3 cos 2 9)d9 


7r 
a 


(5.8) 


Since the partial of F with respect to X does not vanish when X is equal to 
one, the circle docs not extremize the ellipse for the square curvature measure. 



Appendix C: The interpretation of some simple shapes 


In this appendix wc show that our method correctly interprets a number of 
simple shapes, namely an ellipse, a parallelogram, and a triangle. 

Ellipse 

Suppose die ellipse is given by equations (A.2) and (A.3). It is easy to show 
that the area A and perimeter P arc given by 

A — nab 

r2w 

P — (a 2 sin 2 9 + b 2 cos 2 9)%d9 

J 0 

Maximizing A/P 2 is equivalent to minimizing P/\JA. We set b - 
the eccentricity of the ellipse,) and define P/s/A — /(A). We find 

/(\) = — / (—— + x cos 2 e)Ue (c.2) 

yfn Jo X 

By the same argument presented in appendix B, the ellipse will be interpreted 
as a circle provided X = 1 is a minimum of /(X). Changing the variable of 
the integral to <j> = 9 -j- § we find that 

/M = /(i) (C.3) 

which implies that X = 1 extremizes /(X). It is also clear that extrema occur 
in pairs of the form X,l/X. Furthermore, both are stationary points or one 
is a maximum and the other a minimum. Observe that /(X) tends to infinity 
as X tends either to zero or infinity. It follows that a sufficient condition for 
X = 1 to be a global minimum is that all pairs of extrema be stationary points. 
Suppose that this is not the case, and let X 0 be the smallest extremum that is 
not a stationary point. Since /(X) tends to infinity as X tends to zero, Xq must 
be a minimum and 1/Xo a maximum. But 1/Xo is tire largest non-stationary 
extremum, and so, by the same argument as above, it must be a minimum. 
This contradiction establishes the result. 

Parallelogram and triangle 

In section 5 we showed that a skew symmetry is always interpreted as an 
oriented symmetry by our method. In particular, a parallelogram is interpreted 
as a rectangle. By the same argument, a rectangle is a skewed symmetry of a 
square. Hence our method interprets a parallelogram as a square. 

Similar reasoning shows that a triangle is interpreted as a skewed isosceles 
triangle, which is interpreted as a skewed equilateral triangle. The axes of the 
skewed symmetry join a vertex to the midpoint of the opposite side. Hence 
our method interprets a triangle as an oriented equilateral triangle. 


(C.l) 
Xa, (X is 
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Appendix D: The slant and tilt of a skewed symmetry 


In this appendix wc calculate the slant o and the tilt r that correspond to the 
oriented real symmetry that is the interpretation of a skewed symmetry whose 
skew angle is 6. 

As a simple, though instructive, example, consider a rhombus of side a and 
included angle 7 (Figure 8 ). To find the extremizing tilt, wc substitute the data 
from Figure 8 into Eq. 4.24, and find 


sin r cos r 


1 — sin 2 a sin 2 r 


We can rewrite this in the form 


sin(r — 7 )cos(r — 7 ) 
| 2 jl — sin 2 cr sin 2 (r — 7 ) 


P-1) 



|cos 2 (r — 7 ) -f- cos 2 r 

+ (1 — cos 2 a)( 1 — cos 2 (r — 7))(1 — cos 2 r) — 1 


= 0 

P-2) 

We assume first that the first factor is zero. It follows that 


and so 


cos(r — 7 ) = i cos t 


P-3) 


7 — 2 r = MT. P-4) 

Since 0 < r < it and 0 < 7 < tt, there are two possible solutions, namely 


2 

7+ 7T 
2 


P-5) 

P-6) 


Observe that the tilt direction is one of the axes of symmetry o.f the rhombus 
shown in Figure 8 . 

Having solved for the tilt, we now solve for the extremizing slant a using 
Eq. 4.23. Recalling from Eq. D.4 that sin 2 (r — 7 ) = sin 2 r, we find upon 
substitution into Eq. 4.23 




sin 2 a sin 2 r 


1 * 

2 

1 -2 • 2 ) 

} = 

= cos r- 

1 — sin a sin rj 


P-7) 


r-N, 



which we solve to get 


2 cos 2 r 
cos a — - 

sin 2 r 


P- 8) 


The requirement that |cosoj < 1 picks out either Equation (D,5) or (D.6), 
so wc get a unique solution. In this case the skew angle is given by 


« = £- r P-9) 

To summarize, if 5 < 0 we get 

T 

cos a 

and if <5 > 0, 


77 6 

4 ~ 2 
1 -f sin 6 

cos 6 ’ 


(D.10) 


cos a — 


37T 6 

T ~ 2 

cos 6 

1 + sin 6 ’ 


(D . 11 ) 


Tliesc formulae were derived for an equal-sided parallelogram but they will 
clearly apply to the more general case and a rotation through the angles given 
by Equations (D.IO) and (D.ll) will unskew any symmetry. It should be noted 
that r is taken to be zero on the axis of symmetry, as in Figure 8. 

To conclude this Appendix, we consider the case that the second factor in 
Eq. (D.2) is zero and the first factor non-zero. Wc introduce the angle ip by 
analogy with Eq. (2.1) defining orthographic projection: 


tan 2 ip — cos 2 a tan 2 (7 — r) 

Then 


P-12) 


sin 2 (7 — r) — 


cos 2 ^ — r) = 


tan 2 ip 

cos 2 a + tan 2 ip 
cos 2 a 

cos 2 a + tan 2 ip 


(DAO) 


where, without loss of generality, we suppose 0 < ip < 7r/2. Now by 
assumption, the second factor in Eq. (D.2) is zero: 


cos 2 r -f- cos 2 (7 — t) -f- sin 2 <rsin 2 rsin 2 (7 — r) — 1 = 0 P-14) 
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Figure 8 

A typical skew symmetric figure, namely a rhombus of side a and included angle 7 , with 
skew angle 6 = § — 7 . The tilt direction is aligned with one of the axes of symmetry 
of the rhombus, determined by the angle <5. 


Using Eq. (D.13), we find 


/’"‘S 


from which we deduce 


tan 2 rtan 2 xp — 1, 


sin t — cos xp 

cost = n sin ip, M = il, 

since the ranges of the variables can be assumed to be 0 < a < 
t < 7 r ;0 < ip < 7 r/ 2;0 < 7 < 7 r. From Eq. (D.13) we find that 


(B.15) 
7r/2;0 < 


v tan xp cos a 
cos 2 a -f- tan 2 xp ’ 
and, since sin r cos r has the same sign as sin (7 — r) 
Eq. (D.l) we deduce 


sin(7 — r) cos(7 — r) 


v — ± 1 , 

cos(7 — r), \x — 


(D.16) 
v. From 


sin r COST 


1 — sin a sin r 


2 J 


— sin(7 — t) cos(7 — r) 

1 — sin 2 a sin 2 (7 — r ) 2 
v sin xp cos xp 

1 — sin 2 a cos 2 xp^ 


[DAI) 





So far in this Appendix we have only used one of die constraints derived in 
Section 4, namely Eq. (4.24) (from which we derived Eq. (D.l)). We now use 
the second constraint Eq. (4.23), which we can write in the form 

— (1 — sin 2 a 2r)i 

2cos 2 (7 — r) (12.18) 

(1 — sin 2 o sin 2 (7 — r))5 

— (1 — sin 2 a sin 2 (7 — r))^ = 0 

After some algebraic manipulation, we deduce 

cos a — — tan 2 ip, 

from which it follows that cos a is negative, which is impossible in the range 
under consideration. 


2 cos 2 T 

(1 — sin 2 o sin 2 r)£ 
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